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Abstract. We estimate, using high-temperature series expansions, the transition 
temperatures of the spin ^, 1 and | Heisenberg ferromagnet and antiferromagnet in 3- 
dimensions. The manner in which the difference between Curie and Neel temperatures 
vanishes with increasing spin quantum number is investigated. 
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It is well known that in classical spin models, such as as the Ising or classical 
Heisenberg models, on bipartite lattices the critical temperature (if it exists) is the same 
for ferromagnetic exchange (Curie temperature) as for antiferromagnetic exchange (Neel 
temperature). This is a direct consequence of the free energy being an even function 
of the exchange parameter J. It has also been known for some time, but perhaps less 
widely, that for the quantum spin-| Heisenberg model the Curie and Neel temperatures 
are unequal. Early work [TJ |2] put the Neel temperature some 10% above the Curie 
temperature for spin-i, for both the simple cubic (SC) and body-centred cubic (BCC) 
lattices, with the difference decreasing rapidly with increasing S. However these results 
were based on rather short series (six terms) and the critical point estimates contained 
large uncertainties. 

We have re-investigated this question, using substantially longer series (14th order 
for S = ^, 12th order for S = 1, 9th order for S = This is made possible not only 
by the massive increase in computing power now available, but also by the development 
of efficient linked-cluster expansion methods. The reader is referred to a recent review 
jS] for further details of this method. 

The Hamiltonian is written in the form 

H = -jEs. ■ s, -hJ2st - KY.v,st (1) 

where the Sj are spin-S" operators, h and hg are uniform and staggered fields, with 
rji = ±1 on respective sublattices, and the interaction is taken between nearest 
neighbours (ij). J > (< 0) corresponds to the ferromagnet (antiferromagnet). While 
contains the form of the exchange energy for a real spin-S* system, for comparison 
between different S'- values and, in particular, for passage to the classical limit S ^ oo, 
it is convenient to write J = JS{S + 1) and to express critical temperatures in units of 
J/kB. 

The critical temperature kBTc/J is most reliably obtained from the strongly 
divergent 'ordering' susceptibility in zero field: the uniform susceptibility x fo^^ 
the ferromagnet or the staggered susceptibility Xs for the antiferromagnet. High- 
temperature series for these quantities can be derived in the form 

oo 

X,Xs = T.^rK^ (2) 

T- = 

where K = \ J\/kBT and the are numerical coefficients. The uniform susceptibility 
for the spin-| case is known through order K^^ jl], for both the SC and BCC lattices. 
In the present paper we give the staggered susceptibility series to the same order. This 
represents an addition of six new terms to the previously known series jH]. At the same 
time we compute uniform and staggered susceptibilities for the S = 1 case, through 
order K^"^ and K^^ respectively, for both lattices, extending the previous series by five 
terms. We have also calculated the corresponding series for 5* = | through order K'^ . 
For the classical S = oo model the susceptibility series is known through K"^^ [H] and we 
will use this series in our comparison. 
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Tables m 121 and El list the series coefficients, in integer format, for both the SC and 
BCC lattices. The coefficients are positive and appear to be quite regular, suggesting 
that the radius of convergence is determined by the physical singularity on the positive 
real axis (we will return to this point later!) Closer inspection, however, reveals 
some oscillation, reflecting interference from non-physical singularities near the circle 
of convergence. Although we do not base our analysis on this, it is instructive to see 
a ratio plotjTj. We show such a plot in Figure 1. Looking at the SC lattice ffist, it 
is evident that the S = ^ series, in particular, shows a strong 4-term oscillation. This 
results from a pair of singularities on, or near, the imaginary axis, near the circle of 
convergence. The 5 = 1 series are much more regular and, qualitatively, look quite 
similar to the S = oo case. The BCC series are rather regular, even for S* = |. There 
is a 2-term oscillation in all series, which is characteristic of bipartite lattices. Apart 
from the S = ^ (SC) case, resonable estimates of the Curie and Neel temperatures can 
be made visually. Unless something totally unexpected were to occur at higher orders, 
it seems clear that the Neel temperature exceeds the Curie temperature for both S = ^ 
and S = 1 (remembering that the intercept on the ordinate axis is ksTc/J). The very 
similar limiting slope of the different plots is consistent with the universality expectation 
that all quantities diverge with the same exponent 7. 

To obtain more accurate estimates of the critical parameters we turn to Fade 
approximantsj!^. Tables EJ El give estimates of the critical point Kc and exponent 
7, assuming a normal power-law singularity 



obtained from high-order Fade approximants to the logarithmic derivative series. 
Different approximants give quite consistent results and we summarize the overall 
estimates of the critical temperature in Table El The exponent estimates from the 
highest approximants are around 1.42 {S = 1.41 {S = 1). Early studies of the 
S = 00 series also gave values in this range, although the recent long series give lower 
values, approaching the field theory prediction 7 ~ 1.39. Our results are consistent with 
the universality expectation. 

Figure 2 shows plots of critical temperatures ksTc/S^S + 1)J versus 1/S{S + 1). 
The plots appear linear, particularly if the S = ^ points are excluded, and indicate that, 
to a very good approximation 





(3) 




(4) 



where a, b are constants independent of S. Their values are 



SC 



BCC 



a 



b - 



X Xs 

1.4429 
0.288 -0.150 



0.320 



X 



2.0542 



0.174 



Xs 



This linear relation may then be used to obtain reliable estimates of Curie and Neel 
temperature for values 5* > |. 
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With our longer series we are also able, for the first time, to estimate values for the 
amplitudes Cq of the leading singular term Q. This is done in two ways. The first is 
to use the estimates of Kc, 7 obtained previously, form the series for 

{l~K/Kcrx-^Co + ---; K^Kc- (5) 

and evaluate Pade approximants to this series at Kq. The second is to compute the 
series for 

X'^' = C',^\l - K/Kc)-' (6) 

Pade approximants to this series should have a simple pole at Kc with residue KqCq^'^ . 
Both methods give consistent results. We give in Table El our best estimates and error. 
As usual with series analysis, these are not true statistical errors but only confidence 
limits based on the spread of results. As can be seen from Table IHl these amplitudes are 
all of order 1 and show a decrease of some 30% on going from S* = ^ to 5* = 00, with 
the antiferromagnetic amplitude some 5% smaller than the ferromagnetic one. 

The conclusion that the Curie temperature Tq is lower than the Neel temperature 
T/v has a puzzling consequence, as has been remarked on before [2]. Assuming that 
the ferromagnetic susceptibility xi^) ^'Iso has a weak, energy like singularity at —Kn 
[Kn < Kc), as is known to be the case for the Ising model, means that the radius 
of convergence of the series is |-ft"Ar|. Hence the series coefficients must, at some point, 
begin to alternate in sign. To check this point further we follow the procedure of Baker 
et aL[8\, in seeking evidence for a singularity at —Kn in the uniform susceptibility, and 
at —Kc in the staggered susceptibility. To this end we form the series for 



d / d , , 7 

and 



F(K) = —(—\n x(K) ^ ) (7) 

^ ' dK^dK ' Kc-K^ ^ ' 



The first step substracts out the dominant physical singularity from the logarithmic 
derivative series. This series is expected to have a weak singularity at the corresponding 
Neel or Curie point. The final differentiation is to strengthen this singularity In Table 
[7|we show estimates of the location of this secondary singularity and the corresponding 
residue for the S = \ series on the BCC lattice. As is clear, the series F{K) shows 
a consistent pole at if ~ —0.72, consistent with the direct estimate of i^^Ar (Table E}. 
Similarly the series Fs{K) shows a consistent pole at i^' ^ —0.799, consistent with the 
direct estimate of Kc (Table E}. These numerical estimates will, of course, depend on 
the choice made for Kc, K^, 7 in Eqs. ([7j) and (jH)), but are found to be relatively 
insensitive to this choice. We have not repeated this analysis for the SC case or for 
^=1, |. 
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Table 1. Series for x and Xs for spin-|. To avoid fractions a multiplier 4"+^n! for x, 

or 4" + '"(// + 1)! for \ ., lias Ixx^ii iisod. wlicrc^ ii is (lio power of K. 



n 


X 


Xs 




Simple Cubic Lattice S = ^ 





1 


1 
1 


1 





1 


2 


48 


168 


3 


528 


2880 


4 


7920 


59376 


5 


149856 


1478592 


6 


3169248 


42357024 


7 


77046528 


1353271296 


8 


2231209728 


48089027328 


9 


71938507776 


1908863705088 


10 


2446325534208 


83357870602752 


11 


92886269386752 


3926123179720704 


12 


3995799894239232 


198436560561973248 


13 


180512165153832960 


10823888709015846912 


14 


8443006907441565696 


635114442481347244032 




ijody Oentred Oubic 


ijclXXlCe D — 2 





± 


1 
i 


1 


c 
o 


iO 


2 


96 


320 


3 


1664 


8192 


4 


36800 


248768 


5 


1008768 


8919296 


6 


32626560 


367854720 


7 


1221399040 


17216475136 


8 


51734584320 


899434884096 


9 


2459086364672 


51925815320576 


10 


129082499311616 


3280345760086016 


11 


7432690738003968 


225270705859919872 


12 


464885622793134080 


16704037174526894080 


13 


31456185663820136448 


1330557135528577925120 


14 


2284815238218471260160 1 13282648639921512955904 
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Table 2. Series for x and Xs for spin-1. To avoid fractions a multiplier 3""'"^n!/2 
+ l)!/2) has bcH'u used for \ (\,,.) scries, wlicrc n is the power of if. 



Simple Cubic Lattice S = 1 






1 


1 


1 


12 


24 


2 


222 


702 


3 


5904 


26280 


4 


201870 


1184526 


5 


8556912 


63357984 


6 


426905802 


3887604666 


7 


24674144724 


270348199128 


8 


1616505223518 


20988390679758 


9 


118701556096392 


1802403961243776 


10 


9628527879611262 


169418364565523958 


11 


856813238084411136 


17314303199655636792 


12 


82856991914713902402 






Body Centred Cubic 


Lattice D = 1 





1 


1 


1 


16 


32 


2 


424 


1320 


3 


16512 


71136 


4 


819240 


4588968 


5 


50363136 


351263232 


6 


3652143480 


30873601080 


7 


307454670000 


3082065903648 


8 


29310549057000 


343320789071016 


9 


3133368921937824 


42320100429654912 


10 


370060173560963304 


5709664512091086984 


11 


47968071364509850944 


837942419330764322976 


12 


G7.')Gr,427G72.",20r)92:-')4840 
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Table 3. Series for x and Xs for spin-|. To avoid fractions a multiplier 2"+^n!/5 

(2"^ '(;) + I)!/")) lias ])cm used for \ (\ ,.) s(-rics, where n is the power of K. 



Simple Cubic Lattice S = ^ 






1 


2 


1 


60 


60 


2 


1440 


2220 


3 


50136 


106032 


4 


2241660 


6103230 


5 


124125372 


417121164 


6 


8102868414 


32715943017 


7 


613292153184 


2911926450048 


8 


52599376466556 


289263779556198 


9 


5056198898505288 


31792485934519488 




Body Centred Cubic Lattice = | 





1 


2 


1 


80 


80 


2 


2720 


4160 


3 


136448 


283776 


4 


8751600 


23240440 


5 


696028496 


2263139152 


6 


65331028472 


253095247076 


7 


7121212898544 


32175304799424 


8 


879298191968624 


4563926306507096 


9 


121768840349153216 


716734730963510496 



Table 4. Estimates of critical point Kc and exponent 7 (in brackets) from [N, D] 
Fade approximants to the spin-i uniform/staggered susceptibility series. Defective 
PA's are denoted *. 



Spin 5 = i 




Simple 


Cubic 


Body Centred Cubic 


[N,D] 


F(X) 


AF ixs) 


F (x) 


AF ixs) 


[6,7] 


1.1900 


1.0577 


0.7935 


0.7266 




(1.414) 


(1.426) 


(1.416) 


(1.436) 


[7,6] 


1.1925 


1.0611 


0.7935 


0.7266 




(1.432) 


(1.455) 


(1.416) 


(1.435) 


[5,7] 


1.1914 


1.0598 


0.7937 


0.7266 




(1.421) 


(1.440) 


(1.419) 


(1.434) 


[6,6] 


1.1914 


1.0597 


0.7936 


0.7264 




(1.421) 


(1.439) 


(1.417) 


(1.431) 


[7,5] 


1.1931 


* 


0.7939 


0.7267 




(1.438) 




(1.423) 


(1.436) 


[5,6] 


1.1910 


1.0592 


0.7937 


0.7264 




(1.418) 


(1.434) 


(1.418) 


(1.432) 


[6,5] 


1.1901 


1.0583 


0.7936 


0.7264 




(1.411) 


(1.425) 


(1.418) 


(1.432) 
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Table 5. Estimates of critical point Kc and exponent 7 (in brackets) from [N, D] 
Pade approximants to the Spin-1 uniform/staggered susceptibility series. Defective 
PA's are denoted *. 



Spin S = 1 




Simple 


Cubic 


Body Centred Cubic 




^ \x) 




F(X) 


AF ixs) 


[5,6] 


0.38478 




0.26400 






(1.409) 




(1.404) 




[6,5] 


0.38478 




0.26398 






(1.409) 




(1.403) 




[4,6] 


0.38478 


* 


0.26397 


0.25431 




(1.409) 




(1.403) 


(1.405) 


[5,5] 


0.38475 


* 


0.26389 


0.25431 




(1.408) 




(1.398) 


(1.401) 


[6,4] 


0.38467 


0.36541 


* 


0.25410 




(1.406) 


(1.409) 




(1.395) 


[4,5] 


0.38487 


0.36566 


* 


0.25410 




(1.411) 


(1.417) 




(1.395) 


[5,4] 


0.38483 


0.36565 


* 


0.25396 




(1.410) 


(1.417) 




(1.390) 



Table 6. Estimates of the critical temperatures and leading susceptibility amplitudes, 
from Pade approximant analysis. 









SC Li 


\\\kv 










S 


1 

~ 2 


S = 


1 


S = 


3 
2 


S = oo 




F(X) 


AF ixs) 


F(X) 


AF ixs) 


F(X) 


AF ixs) 


X 


Kc 
kBTjJ 
ksTc/J 

Co 


1.192(2) 

0.839(1) 
1.119(2) 
1.26(2) 


1.059(2) 

0.944(2) 
1.259(2) 
1.20(3) 


0.38478(15) 

2.599(1) 
1.2994(5) 
1.11(2) 


0.3656(2) 

2.735(1) 
1.3676(7) 
1.07(4) 


0.195(3) 

5.13(8) 
1.37(2) 


0.190(1) 

5.26(3) 
1.404(7) 


0.69304 

1.4429 
0.9030 


BCC Lattice 




S 


_ 1 

2 


S = 


1 


S = 


3 

2 


S = oo 




F (x) 


AF ixs) 


F (x) 


AF ixs) 


Fix) 


AF ixs) 


X 


Kc 
ksTc/J 
ksTc/J 

Co 


0.7935(3) 
1.2602(5) 
1.6803(6) 
1.15(2) 


0.7266(2) 
1.376(4) 

1.8350(5) 
1.10(3) 


0.2640(2) 
3.788(2) 
1.894(1) 
0.98(1) 


0.2542(2) 
3.934(3) 
1.967(1) 
0.94(1) 


0.1354(10) 
7.39(5) 
1.97(1) 


0.1327(4) 
7.54(2) 
2.009(6) 


0.48680 

2.0542 
0.794 
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Table 7. Estimates of the secondary singularity, a.t —Kn for the uniform susceptibihty 
and at —Kc for the staggered susceptibihty, for the S = ^ models on the BCC lattice. 



[N, D] 


F(K) 










Kc = 0.7936 


7 = 1.416 


Kc = 0.7266 7 


= 1.435 




— -ft'Ar(est.) 


residue 


-i^c(est.) 


residue 


[5,7] 


-0.7160 


0.250 


-0.7993 


0.262 


[6,6] 


-0.7159 


0.250 


-0.7985 


0.260 


[7,5] 


-0.7392 


0.321 


-0.8069 


0.284 


[5,6] 


-0.7189 


0.254 


-0.7987 


0.261 


[6,5] 


* 




-0.7988 


0.261 


[4,6] 


-0.7111 


0.242 


-0.7994 


0.262 


[5,5] 


-0.7107 


0.241 


-0.7985 


0.266 


[6,4] 


-0.7244 


0.275 


-0.7911 


0.241 




1.2 - , , \ / X (S=l/3) 



I I I I I I I I I I I I I I I I I u 

0.1 0.2 0.3 

1/n 

Figure 1. Ratio plots for the uniform 
BCC lattices (as indicated) for 5 = |, 
the variables K = JS{S + 1)/A:bT, or e 
are the coefficients of the iiT-series 



2.2 



1.8 



1 1 1 1 1 
. BCC 






^ - 


















0.1 0.2 0.3 



l/n 

and staggered susceptibilities for the SC and 
1, oo. The ratios are defined for the series in 
quivalently r„ = a„/(S'(S' + l)a„_i) where a„ 
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